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Abstract. In this paper, we investigate the strong gravitational lensing by the stationary, 
axially-symmetric black hole in Kerr-Taub-NUT spacetime in the strong field limit. The 
deflection angle of light ray and other strong deflection limit coefficients are obtained numer- 
ically and they are found to be closely dependent on the NUT charge n and spin a. The 
magnification and the positions of the relativistic images are computed. The caustics are 
studied and the results show that these caustics drift away from the optical axis, which is 
quite different from the Schwarzschild black hole case. Moreover, the intersections of the 
critical curves on the equatorial plane are obtained and it is shown that they increase with 
the NUT charge. These results show that there is a significant effect of the NUT charge on 
the strong gravitational lensing. 
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1 Introduction 

Deflection of light ray in a gravitational field is one of the consequences of Einstein's general 
theory of relativity. The phenomenon is referred to as gravitational lensing and the object 
causing a detectable deflection is called a gravitational lens. Researches on the gravitational 
lensing show that the value of the deflection angle of light depends on the observer-lens 
distance and the nature of the lens. This result implies that the gravitational lensing can 
provide us with the information about the distant stars. It also can help us to test the 
exotic objects in the universe and to estimate the values of the cosmological parameters [1]. 
Moreover, it could provide a profound verification of alternative theories of gravity in the 
strong field regimes [2-5], and an effective detection of gravitational waves [6, 7] and the 
cosmic censorship hypothesis [8, 9]. 

Gravitational lensing in weak field approximation has been developed in [10-12]. In the 
weak field approximation, one can study the properties of the ordinary stars and galaxies. 
However, when a compact object (like a black hole) plays the role of the lens, there will be 
notable phenomenon very near the object, which can be described by the strong field limit 
or the strong deflection limit rather than the weak field approximation. The studies of the 
strong field limit lensing due to different black holes have received considerable attention 
in recent years, which suggests that we can extract the information of the black hole from 
lensing. Studies on it can be traced back to the work [13, 14] of Darwin, as well as to that 
of Frittelli, Kling and Newman [15] and that of Virbhadra and Ellis [16], where the authors 
presented a definition of an exact lens equation without reference to black hole background 
spacetime, and then the exact lens equation was constructed for the Schwarzschild black 
hole spacetime. It is worth noting that the Virbhadra-Ellis lens equation is effective for 
the case that the observer and the light source must be far away from the lens so that the 
gravitational fields there can be described by flat metric. Later, the work was extended to the 
Reissner-Nordstrom black hole lensing [17] and the naked singularities lensing [8, 9], which 
showed that the strong gravitational lensing by naked singularity is very different from that 
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of black hole. With the lightlike geodesic equation, lensing in a spherically symmetric and 
static spacetime was considered in [18]. 

Based on the Virbhadra-Ellis lens equation, Bozza et al, [19] proposed an analytical 
method for obtaining the deflection angle in the strong gravitational field case and the result 
shows that the deflection angle diverges logarithmically as the light rays get close to the 
photon sphere of a Schwarzschild black hole. In [20], Bozza proved that the result is also held 
for other static spherically symmetric black hole lensing. He still extended the methods to 
study the spinning black holes [21-23]. Moreover, Eiroa et al, have studied the gravitational 
lensing by the Reissner-Nordstrom black hole and the braneworld black hole [24-26] . Strong 
gravitational lensing by other black holes and wormholes were extensively investigated in 



On the other hand, the Kerr-Taub-NUT (KTN) black hole [48, 49] is a remarkable 
solution of Einstein-Maxwell equations for electro-vacuum spacetime possessing with grav- 
itomagnetic monopole. The KTN black hole carries three parameters, the mass M, the 
spinning parameter a and the NUT charge n known as the "gravitomagnetic mass". The 
presence of the NUT charge gives rise to a black hole solution with some fascinating proper- 
ties. For example, there exists no curvature singularity but conical singularities on its axis 
of symmetry, which is caused by the gravitomagnetic analogue of Dirac's string quantization 
condition [50]. Furthermore, the influences of the NUT charge or gravitomagnetic mass n 
on the spacetime structure, thermodynamics and geodesic equations, as well as the particles 
collision were studied in [51-59]. 

The gravitomagnetic lensing in nonrotating NUT spacetime has been studied in [60]. 
The results show that the NUT charge influences the gravitational lensing through the null 
geodesies. The possibility to detect the gravitomagnetic masses with the next generation of 
microlensing experiments was analyzed in [57, 61, 62]. These results suggest that the influence 
of the NUT parameter will be stronger in the vicinity of compact gravitating objects with 
small radius. So the purpose of this paper is to study the effects of the NUT charge n and 
spin a on the strong gravitational lensing by the rotating supermassive KTN black hole. 

The paper is structured as follows. The second section is devoted to the derivation of 
the first order differential system for the geodesies in the background of the KTN black hole. 
In Sec. 3, we study the lensing equation in the equatorial plane and calculate numerically the 
strong field limit coefficients and the deflection angle. And the differences of these coefficients 
between the Kerr black hole and the KTN black hole are analyzed. In Sec. 4, we study the 
quasi-equatorial lensing by the KTN black hole, where the precession of the orbit for small 
declinations, and the magnification of the images are obtained. The critical curves and caustic 
structure are analyzed in Sec. 5. And a brief discussion is given in Sec. 6. 

2 Geodesies in Kerr-Taub-NUT spacetime 

The general KTN black hole solution of the Einstein field equations is described by the metric 



[27-47]. 



[48, 49] 




(2.1) 
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where 



p2 = ^2 _^ (n + acosef, (2.2) 
A = r"^ -2Mr + a'^ -n^. (2.3) 

Here M is the gravitoelectric mass of the black hole and n is the NUT charge (or gravito- 
magnetic mass). The angular coordinates 9 G (0, vr) and (j) £ (0,27r). In the particular case 
n = 0, the solution (2.1) coincides with the Kerr black hole solution. And in the case a = 
and n = 0, it recovers the Schwarzschild black hole solution. 

The event horizon of the KTN black hole is located at the biggest root of the equation 
A = 0, i.e., 

Th = M + \/m2 + n2 -a2. (2.4) 

For the case > + n^, there will be no event horizon and a naked singularity will appear, 
which spoils the causality of the spacetime and is forbidden according to the Penrose's cosmic 
censorship conjecture. So, we here only consider the case < + n^. Another important 
surface of the black hole is the ergosphere, which corresponds to gu = 0. The outer one is 

res = M + \/m2 + n2 - a2 cos^ 6. (2.5) 

Note that we always have res ^ fh and the inequality is saturated at = 0, vr. 

The geodesic equations can be solved with the Hamilton-Jacobi method for the metric 
(2.1) and the first-order differential system for the geodesic is given by (the motion of a 
charged particle around the KTN black hole immersed in an external magnetic field can be 
found in [56]) 

= ±Vr, (2.6) 
= ±Ve, (2.7) 



E 



sin 9{a + 2ncsc of -An{n + a{l + cos 9)) 



(2n cos 9 csc^ 9 — a)L 



+ 



(r2 + a? + v?)[E{r'^ + a? + n?) - aL] 
A ' 

(a - 2n cot 9 esc 9)E + L csc^ 9 + 



a[£'(r2 + + n'^) - aL] 
A ' 



(2.8) 
(2.9) 



where the dot denotes the derivative with respect to the affine parameter and R, are 



R = {E{r'^ + a^) - aL)2 + n^E{2a'^E - 2aL + Ein"^ + 2r'^)) 

- [{L - aEf + m2(n2 + r^) + K] A, 
e = /C - cos^ 9[a'^{m^ - E^) + csc^ 9] 

+2n cos9[2aE^ - am^ + 2EL csc^ 9] - An^E"^ cot^ 9. 



(2.10) 



(2.11) 



Equations (2.6)-(2.9) are the first-order geodesic equations for a massive particle m? = 1 and 
for a photon rn? = 0, respectively. /C is a separation constant of motion. The constants E 
and L are the conservation of energy and orbital angular momentum per unit mass of the 
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motion and they correspond to the Kihing fields dt and 5^, respectively. For the case n = 0, 
Eqs. (2.6)-(2.9) just describe the motion of a particle in the Kerr black hole spacetime. 
Moreover, we can express the lightlike geodesies in the following form 

dr dO 

^7=, (2.12) 



, [^(r^ + + n^) - aLl , /"^ L csc^ 6* - (a - 2n cot 61 esc 6l)£; ,^ 
A(p = a — ^ -dr+ ^ — ^—d9, 2.13 

^ r (r^ + + n^)[E{r^ + + n^) - aL] 

(a - 2n cos 9 csc^ 9)L - ^'[sin^ 9{a + 2n csc^ 9)"^ - An{n + a(l + cos 9))] 

d9. 



(2.14) 

Note that the NUT charge n indeed influences on the null geodesies, which will be shown to 
have significant effect on the black hole lensing. 

As we know, in the real world, the graviational field far away from a compact object 
or even a black hole is very weak and can be described by a fiat metric. So the compact 
object only significantly infiuences the motion of particle in the neighborhood of the object. 
On the other hand, we take the small value of n/2M for the quasi-equatorial approximation 
with 9 ^ ^. Then the light ray trajectory can be regarded as a straight line at infinity both 
in equatorial plane and quasi-equatorial approximation. From this view, one can identify 
the approximate light ray with three parameters ■0o) u and h. The first one ■0o represents 
the inclined angle that the incoming light ray forms with the equatorial plane. The second 
one u is an impact parameter of the projection of the light ray trajectory in the equatorial 
plane. And the last one h describes the height between the point of the projection closer to 
the black hole and the trajectory. Then following [21, 30, 33], if the observer is located at 
{fo, "^o) in the Boyer-Lindquist system, one can define two celestial coordinates Ci and C,2 
for an image. The coordinate Ci denotes the observable distance of the image with respect 
to the symmetry axis in direction normal to the ray of sight and the coordinate C2 measures 
the observable distance from the image to the source projection in the equatorial plane in 
the direction orthogonal to the ray of sight. With the help of Eqs. (2.6)-(2.9), it is easy to 
express the two coordinates C,i and C2 in the following form: 



dr 



^d^ 

C2 = ro-r 



r, To^oo 



dr 



(L + 2ncosi?o)sin~^t?o, (2.15) 
hsm^o, (2.16) 



r, To— >oo 



where we have taken the choice of S = 1. Considering the quasi-equatorial case (i.e., "!?(, = 
7r/2 — tpo with small tpo) and Ci = we can obtain the orbital angular momentum L and 
Carter constant /C, in terms of u and ipo'- 

L = ucosipo — 2nsimpo, (2-17) 
K. = cos^ ipo + (u^ — g^) sin^ ipo — 4nsin ^^0(0 — ncos + nsinipo)- (2.18) 

It is clear that, when n = 0, the result will reduce to the Kerr black hole case. 
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3 Equatorial black hole lensing 



We devote this section to study the equatorial lensing by the rotating KTN black hole. The 
effects of the NUT charge n and spin a on it will also be investigated. 

3.1 Deflection angle 

In this subsection, we would like to consider the strong field lensing by the KTN black hole 
for the case that both the observer and the source lie in the equatorial plane {9 = tt/2) of the 
KTN black hole and the whole trajectory of the photon is also limited on this plane. Here, 
we can adimensionalize the metric (2.1) in terms of the Schwarzschild radii 2M by defining 

t t/2M, r r/2M, a a/2M, n n/2M. (3.1) 

Then, the metric is reduced to 

ds^ = -A{r)dt^ + B{r)dr'^ + C{r)d^'^ - D{r)dtd()), (3.2) 

with the metric coefficients given by 

2 2 

= 5 3.3 

- (3.«) 

In the equatorial plane, we can express the first-order geodesic equations (2.6)-(2.9) for the 
photon {m? = 0) in terms of the metric coefficients A{r), B{r), C{r) and D{r), which read 



ACE - 2DL 



(3.7) 



CE^ - DEL - AL^ , , 

* = ^'i—mcTiwr- 

^ = 0, (3.9) 
2DE + 4AL 

For simplicity, we take the choice E = 1. On the other hand, we could rewrite (3.8) as 

r' + K=ff = 0, (3.11) 

where the effective potential reads 

V B{4AC + D^) J 
_ 2aL{2n? + r) -{-n? + r - r"^)!? - a'^{'in? + r + r^) 
= (n2 + r2)2 ^^-^^^ 
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At the minimum distance of photon trajectory, where T4ff = 0, one could get [21] 



L/ — U — : 

2An 



_ -a(ro + 2n^) + (rp + n'^)y/ro{ro - 1) + - 

- ro - 

where we have used Eq. (2.17). The subscript "0" represents that the metric coefficients are 
evaluated at rg. Here, the sign before the square root has been chosen to be positive. Thus, 
a prograde photon orbit is related to a > 0, and a retrograde one is related to a < 0. 

With the expressions (3.8) and (3.10), the deflection angle for the photon coming from 
infinity can be written as 

a{ro) = (P{ro)-T^/2. (3.14) 

Here, denotes the period of the angular coordinate (p. As we know, there exists a conical 
angle on its axis of symmetry. However, for small n, we have ~ 27r for approximation. 
Thus the total azimuthal angle (p{ro) is given by 



oo 



Hr.) = 2l ^^BmD^2LA) ^^^^^^ 



■0 VaACTd^ ^sgn{Ao) {CAo - ACo + L{ADo - AqD)) 

where, sgn(X) denotes the sign of X. Note that the metric coefficient Aq changes its sign 
when ro < res and we have taken this into consideration. With a detailed examination, we 
can see that the deflection angle a(ro) increases with the decreasing of the parameter rg. At 
a certain value of tq, we may have a(ro) = 27r, which means that the light ray will make a 
complete loop around the black hole before reaching the observer. Let ro decrease further, 
the light ray may make more than one complete loop, and when tq approaches the photon 
circle radius Tc (we will explain it later), the deflection angle a(ro) will be unboundedly large 
and the photons will be captured by the black hole. 

Following the method developed by Bozza [20] , we can find the behavior of the deflection 
angle when photons get very close to the radius rc. We first define two new variables y and z 

y = Air), (3.16) 
z = , 3.17) 

1 - yo 

where yo = Aq. With the two new variables, we can express the total azimuthal angle as 

(t){ro)= [ R{z,ro)f{z,ro)dz, (3.18) 

with 







Riz,r.) = '-^y^^B±2£S, (3.19) 

/(z,ro) = ^= ^ =. (3.20) 

'-^^ [CAo - ACo + L{ADo - AqD)] 
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These metric coefficients without the subscript "0" are evaluated at r = — yo)z + yo)- 

Note that the function R{z, tq) is regular for all values of z and ro, while f{z, rg) diverges at 
z = 0. Thus, the integral (3.18) can be separated into two parts 

<P{ro) = ^R{ro) + <PD{ro), (3.21) 

with the divergent part 

Mro)= f R{0,Xc)fo{z,ro)dz, (3.22) 







and the regular part 



(pniro) = / g{z,ro)dz 
Jo 



(3.23) 



with g{z,rQ) = /^(z, ro)/(z, tq) — R{0,rc)fo{z,rQ). In order to find the divergence of the 
integrand, we expand the argument of the square root of /(z, rg) to second order in z and 
the function fo{z,rQ) is 

fo{z, ro) = , \ (3.24) 

y/pz + qz^ + 0{z^) 



where 



p = sgniAo)^^^r^(^AoC'o - A',Co + L{A',Do - ^o^o)) , (3.25) 
q = sgn(Ao) ^^^~,^;]' hcoC'oA'i + (CoC^ - 2C^)AoA'o - CoC'^AoA'^ 



+L [AoCoiAM - A'oD'^) + 2A'oC'oiAoD'o - A'^Dq)] ) . (3.26) 

When the coefficient p vanishes and the leading term of the divergence in /o is z~^, we will 
obtain an unlimited deflection angle. Here let us give some notes on the photon sphere (for 
a = 0) or photon circle (for a 7^ 0). As we know that, for a static spherically symmetric 
spacetime, there exist several definitions of photon sphere [16, 63]. One definition states 
that it is a timelike hypersurface if the Einstein bending angle of a light ray is unlimited 
when the closest distance of approach coincides with the photon sphere. Extending this case 
to the axisymmetric spacetime, there will be a photon circle rather than photon sphere in 
the equatorial plane. When the closest distance of photon orbit coincides with this photon 
circle radius Tc, we will get an unlimited deflection angle. Therefore, solving p = 0, we will 
obtain the radius Tc of the orbit, which is determined by the following equation (for a detailed 
definition of photon sphere in a stationary and axisymmetric spacetime, we refer readers to 
[47]) 

A,C', - A',Cc + L,{A',D, - A,D',) = 0. (3.27) 

Here, the prime indicates the derivative with respect to r and the subscript "c" represents 
that the metric coefficients are evaluated at r = re- In fact, we can obtain the photon circle 
radius Vc through the effective potential Vcs (3.12) by solving Vcs = and = 0. 
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Here for the KTN black hole case, the radius Tc is determined by the following equation 

4r^ - I2rl + 3(3 - - 8(a^ - ^n^)rl 

+ 2n2(18n2 - IGa^ - 3)r^ + An'^{2a^ - ?,n^)rc + = 0. (3.28) 

Obviously, this equation is of 6th degree in the radius r^. So, for the fixed values of a and 
n, we will obtain six values for r^. However, there are always two imaginary values and two 
negatives, which are unphysical. And the two remaining ones are larger than r^. From the 
physical viewpoint, the larger one of the two is the radius of the retrograde circular photon 
orbit with a < and the smaller one is that of the prograde circular photon orbit with o > 0. 
Thus, consider that |a| = a and \a\ = —a together, is computed numerically and is plotted 
in Fig. 1 for different values of n. The black solid line with n = describes the case of the 
Kerr black hole, which has a analytical expression = l + cos(| arccos(ib2|a|)) corresponding 
to the radius of the prograde and retrograde circular photon orbits, respectively. For n = 
and a = 0, we get the radius of photon sphere rpg = | for the Schwarzschild black hole. We 
can see that, for a fixed spin a, the radius increases with the NUT charge n. Note that 
the light ray trajectory is a straight line at infinity in equatorial plane, so the value of n/2M 
could take a large value for this equatorial black hole lensing. It is also clear that, for a 
positive a, photons are allowed to get closer to the black hole, and enter even the ergosphere 
at some values of the spin a. We calculate the critical spin Ccr, where the circular radius Tc 
coincides with the ergosphere in the equatorial plane 

Vl + 4n2 

». = ^^. (3.29) 

It reduces to the Kerr black hole case a^r = for n = 0. The behavior of the critical spin 
Ocr is shown in Fig. 2. The shadow area on the top of Fig. 2 represents the case of naked 
singularity. The remaining area is divided into two regions by the critical curve. Below the 
curve, we have sgn(^c) = 1 and sgn(Ac) = —1 for another region. It is also clear that, from 
Fig. 1, the prograde photons always have smaller values of than the retrograde ones. So, 
we may conclude that the retrograde photons are captured more easily than the prograde 
ones. It is also clear that, for a fixed spin a, the Kerr black hole has the minimum value of 
the radius Tc, which implies that the photons are more easily captured by the KTN black 
hole than by the Kerr black hole. 

For the case ~ Tc, the deflection angle can be expanded in the following form [20] 

u — 

a{u) = -a\og{ l) + h + 0{u-Uc). (3.30) 

The coefficients Mc, a and h are given by 



Uq — L\rQ=rci (3.31) 

(3.32) 



i2(0,rc) / 2AcSe 



2^ \ A,C'l- A'lC + u, {A'lD, - A,D>>) ' 
"" + '« + '^°^(uJ^J('j:^2.,^.))' ^'-''^ 



where Qc = q\r=rc- I^i order to obtain the coefficient b^, we expand (pfji^o) at 

n=0 



CO 

zTn ^- Jo 



dr^ 



dz. (3.34) 

ro=rc 
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n/2M=0.0 

n/2M=0.4 

n/2M=0.8 

n/2M=1.2 

n/2M=1.6 
n/2M=2.0 
n/2M=2.4 



-2 



-1 




a/2M 



Figure 1. The radius Tc as a function of the spin a for the Kerr black hole (black solid line) and the 
KTN black hole. 




nl2M 



Figure 2. The critical spin Ocr vs the NUT charge n. The shadow area denotes the case of naked 
singularity. 



Thus, we get 



bn = 0ij(rc) = / giz,rc)dz, 
Jo 



(3.35) 



which can be obtained numerically. 
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Figure 3. Variation of tlic minimum impact parameter Uc with the spin a of the Kerr black hole 
(black solid line) and the KTN black hole. 
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Figure 4. Strong deflection limit coefficient a as a function of the spin a of the Kerr black hole (black 
solid line) and the KTN black hole. 



Fig. 3 depicts the minimum impact parameter Uc as a function of the spin a for different 
values of the NUT charge n. We can see that Uc has the similar behavior as Tc- We also can 
see that the Kerr black hole has smaller value of Uc for a fixed spin a than that of the KTN 
black hole. The coefficients of the strong deflection limit a and b are illustrated in Figs. 4 
and 5. It is easy to obtain that, for a fixed NUT charge n, a grows with the spin a, while 
b decreases. And both the coefficients diverge at a = ^ \/l + 4n-^ , which corresponds to the 
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Figure 5. Strong deflection limit coefRcient 6 as a function of the spin a of the Kerr black hole (black 
solid line) and the KTN black hole. 
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I a/2M 

-2-10 1 2 

Figure 6. The deflection angle a(u) vs the spin a of the Kerr black hole (black solid line) and the 
KTN black hole for u = Ur + 0.0025. 



extremal black hole. The divergence of the coefficients implies that, in the strong field limit, 
the deflection angle no longer represents a reliable result. 

With a and b, we can obtain the deflection angle a{u). The change of a{u) with spin a 
for u = Uc + 0.0025 is displayed in Fig. 6. From it, we obtain the results: (1) For the Kerr 
black hole (described by the black solid line in Fig. 6), a{u) monotonically increases with 
the spin a, and a{u) of prograde photon is larger than the retrograde one. (2) For the KTN 
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black hole with the NUT charge n 7^ 0, the change of the deflection angle a{u) with the spin 
a has the same behavior as the Kerr black hole. (3) For a fixed value of the spin a less than 
some negative value, the deflection angle a{u) of the Kerr black hole is smaller than that of 
the KTN black hole; and for a spin a greater than this negative value, a{u) of the Kerr black 
hole is larger than that of the KTN black hole. 

3.2 Lens equation 

In order to describe the lensing, several lens equations were introduced (such as [64-66]). 
These equations principally differ from each other for the different choices of the variables. 
In [67], the author gave a detailed comparison of each lens equation and he suggested that 
the Ohanian lens equation is the best approximate lens equation. So, in this subsection, we 
would like to give a brief introduction to the equatorial lens equation. General, the optical 
axis is defined as the line joining the observer and the lens. Setting the black hole in the 
origin, we denote the angle between the direction of the source and the optical axis by 7. 
Then the case 7 = corresponds to that the source, lens and observer are perfectly aligned. 
From the lensing geometry, the angle 7 can be expressed as 

7 = -a(6l) + 61 + ^ mod 27r, (3.36) 

where the impact angle is 

- u ODoL „„x 

ASSIST- 

Although the spacetime has a conical singularity, the lens equation (3.36) is still held, because 
that all azimuth angles are measured in coordinate (p with period = 27r. Dls measures 
the distance between the lens and source, -Dql for the observer and lens, and Dqs for the 
observer and the source [67]. A usual relation between them Dqs = -Dql + -Dls is held. The 
lens equation can be reexpressed as 

7 = ^OL + -DLS g _ ^^^^ 2vr. (3.38) 

Since the angle 7 € [— vr, vr], the source and the observer could be on the opposite side or on 
the same side of the lens. 

Consider that 9 = u/Dql 1, we have 7 ~ —a{9). Moreover, one may find 

0O = -^(l + e„), (3.39) 

where = e(^+7-2n7r)/a ^^^^^ ^ -g number of loops done by the photon around the black 
hole. Expanding the deflection angle a{9) around 9^ to the first order, we get 

da I 
'89 



a{9) = aie^J + ^ {9 - a°) + 0{9 - 9°J' mod 27r 



~ _^ _ ^^^(e - 9^) + 0{9 - 9lf mod 2Ti. (3.40) 

uc 

Neglecting the higher order terms and plugging this result into the equatorial lens equation 
(3.38), we obtain the position of the n-th relativistic image 

^nC^^^fl- '-^-i^^^ + '^^^^ V (3.41) 

From (3.41), it is easy to find that the correction is much smaller than 9^ for Dql ~ -Dls S> 1- 



- 12 - 



4 Quasi-equatorial black hole lensing 

It is known that, in an axisymmetric spacetime, even if the central ray of the hght bundle is 
in the equatorial plane, most of the rays in the bundle will leave the equatorial plane. Thus, 
in this section, we would like to investigate the quasi-equatorial lensing by the KTN black 
hole. It was suggested by Bozza [21] that, in order to study the quasi-equatorial lensing, a 
two-dimensional lens equation is needed. Accordingly, an additional parameter ip = Tr/2 — -d, 
which measures the declination between the orbit plane of the light ray and the equatorial 
plane, should be included in. It is also worth to mention that, we here adopt small n/2M 
approximation and '0 ~ 0, so the technique for the Kerr black hole case is also applicable for 
this case. 



4.1 Precession of the orbit 

The precession of the orbit at small declination is an important ingredient for the study of 
quasi-equatorial lensing. Here, we only consider the small n case. Similar to the Kerr black 
hole, the evolution equation for ip as a function of the azimuthal angle (p for the KTN black 
hole approximatively reads 



^^±u;{4>)\/i;^-i^^, (4.1) 

where 



u = y^u^-a^, (4.2) 

= + (4-3) 

- -n^ + r{r - 1) ,^ ^, 

ojU) = u— TTz. (4.4) 

It is obvious that the parameter Cl'((/)) depends on the NUT charge n and spin a. When n = 0, 
uj{(l)) will reduce to the Kerr black hole case [21]. The solution of (4.1) is 

V'(</') = V'cos((^ + (/)o), (4.5) 

with = Jq Lj{(j)')d(j)' and (j)Q a constant. Since the photon comes from infinity and returns 
to infinity, the deflection angle is 



<pf = 2l u;{r)^dr = J^R^{z,ro)f{z,ro)dz, (4.6) 



where 



RUz,ro)=ioir)R(z,ro). (4.7) 

The quantities R{z,rQ) and /(z,ro) are given by Eqs. (3.19) and (3.20), respectively. Since 
there are no singularities in uj{r), it can be absorbed into the regular function R{z,rQ). 
As a result, we can apply the same technique used in section 3 for this case. Thus, the 
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Figure 7. The coefficient a vs the spin a for different values of the NUT charge n. 



representation of in the strong field limit reads 



1 



4'f = —a In 

. i?c.(0,rc) 
a = . 

h = hn + aln 



with 



bR 



Rujiz,ro)f{z,ro) - RujiO,rc)foiz,ro) 



dz. 



(4.8) 
(4.9) 
(4.10) 

(4.11) 



The quantities a and b are evaluated numerically at the circular radius Vc and their behavior 
is shown in Figs. 7 and 8. From them, we can find that, for the Kerr black hole, we always 
have a = 1 for all spin a. However, for nonvanishing charge n, both quantities a and b are 
found to increase with the spin a. And both of them diverge at a = 1/2. Variation oi (j)f 
with the spin a is presented in Fig. 9 with u = Uc + 0.0025. For fixed charge n, it increases 
with the spin a. And for fixed spin a, it decreases with the charge n. 

4.2 Lensing at small declination 

In this subsection, let us study the positions and magnification of the images at small decli- 
nation. Here we set the source and observer at heights hs and ho, respectively. In order to 
ensure the small declination condition during the whole trajectory of the photon, we assume 
that the following relation is satisfied 



iho,hs) <u< {Dol,Dls)- 



(4.12) 
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Figure 9. Variation of the values of 0/ with the spin a for fixed u = Uc + 0.0025. 



With the small declination hypothesis, the polar lens equation connecting the positions of 
the source and the observer can be expressed as [21]: 

hs = ho - - (pOL + Dj^s)C - , (4.13) 
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where S = sin<j)f and C = cos(pf. The inchnation ■00 is related to the heights of the observer 
and the source. For the n-th image, the inclination ipQ and. the height h^i are 



hr, 



u{hs + hpCn) - hoDi^sSn 
DolDi^sSu - u{DoL + Di^s)Cn ' 
u{hsDoL - hoDLsCn) 



(4.14) 
(4.15) 



where Sn and C„ are the values of S and C evaluated at (f)f = (pf^n- It is clear that in 
the neighborhood of (j)f = kn, the denominators of the ipo^n and hn will vanish. Thus, the 
position of the caustic points is naturally determined by the denominators of the Vo.n and 
hm i.e., 



Kij) = DolDlsSu - uiDoL + l?Ls)C„. 



(4.16) 



So, the caustic points are located at (j)f = kir. Solving it, we could obtain the angular 
positions of the caustic points: 



lk = -b+ -{b- kn). 
a 



(4.17) 



For each k, there exist one caustic point for the prograde photons and one for the retrograde 
photons. The case k = 1 corresponds to the weak field caustic points and others to that in 
the strong field limit approximation. 

Next, we would like to study the magnification of the images near the caustic points. 
The magnification is defined as the ratio of the angular area element of the image and that 
of the source without lens, 



= 

The Jacobian determinant reads 

|J| : 

where 



\J\ 



dhs 



aDoi.\K{^)y 



,(fe+7)/a 



(4.18) 



(4.19) 



(4.20) 



The number of loops n = made by the photon has been coded in 7, so 7 here can take 
an arbitrary value. And he angular position of the source is described by 7 mod 27r. It is 
also worth pointing out that different values of 7 differing by a multiple of 2tt represent the 
same source position with respect to the lens but reached by photons performing a different 
number of loops. 

We show the positions of the first five relativistic caustic points in Fig. 10. When the 
charge n = 0, it describes the result for the Kerr black hole and it is consistent with the 
result in [21]. When n = and a = 0, it recovers the result of the Schwarzschild black hole 
and we have 7^ = —{k — l)7r, which means all the caustic points for the Schwarzschild black 
hole are located at the optical axis. While for the Kerr black hole with or without the NUT 
charge n, all caustic points of the source are not on the optical axis. We can also see that 
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Figure 10. The angular position of tlic first five relativistic caustie points. From below to above 
k = 2,3,4,5,6. 



the caustic points are anticipated for negative spin a and delayed for positive spin a. It is 
obvious that the caustic curves can move very far from the optical axis at large values of 
the spin a. Moreover, we can obtain the result that the NUT charge n has strong impact on 
the prograde photons than the retrograde ones. Prom Eq. (4.19), we find that the Jacobian 
determinant will diverge in the caustic points. Thus, in order to describe the magnification 
of the enhanced images created by a source near the caustic points, we expand (4.16) around 
7fc and retain the first term 

a 

Then the magnification of the enhanced images is expressed as 

^enh ^ POL + I^LS)^ /2.(a) 
' ^k^LS^ 17-7.1' ^ ' 

fik{a) = uuce^^{a)a~^ , (4-23) 

e^Jo) = e(^-*^")/^ (4.24) 

The quantity flk denotes the magnifying power for the KTN black hole close to the caustic 
points. The behaviors of Uk for k = 2,3,4,5 are shown in Fig. 11. From it, we can get 
the results that first decreases linearly with and then grows with the spin a. And each 
Uk shows a divergence for the extremal black hole case. However, we should keep in mind 
that the result that the behavior of the magnifying power around an extremal black hole is 
inaccurate because the strong field limit approximation there breaks down. It is also clear 
that the magnification falls rapidly with the increasing of k. 

It is obvious that the shapes of p.^ shown in Fig. 11, remain more or less the same for 
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Figure 11. Behaviors of the magnifying power /i^ for fc = 2, 3, 4, 5. 



different values of k. And for two consecutive numbers k and /c + 1, we get 

Since d = 1, the Kerr black hole has a constant ratio fik+i/jj'k = e~'^ . The numerical results 
of the ratio for the NUT charge n ^ can be found in Fig. 12. We can find that, for fixed 
charge n ^ 0, the ratio increases with the spin a. The ratio decreases with the charge n for 
fixed spin a with small value, and the result for fixed high spin a is the reverse. 

5 Critical curves and caustic structure 

It is known that, for the Schwarzschild black hole lens, if the source, lens and the observer 
are strictly aligned, then a large Einstein ring and two infinite series of concentric relativistic 
Einstein rings will appear [20]. However, for the rotating black hole lens [21, 33], the result 
shows a big difference. The nontrivial caustic structures appear, and the caustics drift away 
from the optical axis and acquire a finite extension. Especially, for a black hole with high 
spin a, only one image rather than two infinite series of relativistic images will be observed. 
In this section, we would like to study the effect of the NUT charge n on the critical curves 
and the caustic structure. 



(4.25) 
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Figure 12. The ratio fik+i/fJ-k vs. the spin a. The black sohd hne has a constant value e 



For the KTN metric, the intersections of the critical curves with the equatorial plane in 
the quasi-equatorial approximation and in the strong deflection limit are 

V aiJoL-i-'LS / 

where 

e^ = 7^(l + e,J. (5.2) 

Here, we suppose that the gravitational field of the supermassive black hole at the center of 
our Milky Way can be described by the KTN metric. The mass of the supermassive black 
hole is estimated to be M = 2.8 x W^Mq, and the distances are assumed to Dql = 8.5 kpc, 
Dls = 1-0 kpc. Then the numerical results for 6'^ are plotted in Fig. 13. From it, we can 
find that these curves are very close to each other and the critical points are close to the 
optical axis 9 = for positive spin a and farther for negative spin a. 

From above discussion, we get a result that the caustics for a rotating black hole exist 
a nonvanishing extension. Here we would like to study it and draw sketch maps for it. It is 
known that the intersection of the k-th caustic with the equatorial plane is determined by 
7fc(— |a|) and 7fc(|a|). We plot the first six caustics in Fig. 14 and Fig. 15 for fixed spin a and 
fixed NUT charge n, respectively, seen from the direction of the spin. In these figures, the 
centre of the circular sectors stand for the location of the black hole lens. The line connecting 
the centre and O (observer) is the optical axis. The circular arc denotes the angular positions 
7fc of the source. And the distances between the source, lens and observer are ignored in these 
figures. From them, we can find that the nonrelativistic caustic 71 stays close to the optical 
axis. And other relativistic ones drift in the anticlockwise direction for fixed spin a, and in 
the clockwise direction for fixed NUT charge n. In Fig. 14, the relativistic caustics 72, 74 
and 76 are shifted to below the optical axis, and 73, 75 are shifted to above it. While, in 
Fig. 15, these relativistic caustics are shifted in the opposite direction. It is also clear that, 
for fixed NUT charge n, the caustics get larger and farther from their initial position on the 
optical axis when spin a and k increases. However, for fixed spin a, the caustics keep its 
extension and get closer from their initial position. 
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Figure 13. Intersections of the strong deflection limit critical curves with the equatorial plane for 
k = 2,3,4,5. 
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Figure 14. The first six caustics for the KTN black hole with fixed spin a = 0.12, marked by the 
arcs between 7fe(— |a|) and 7fc(|a|). 



6 Discussions and summary 

In this paper, we studied numerically the equatorial and quasi-equatorial lensing by the 
stationary, axially-symmetric KTN black hole in the strong field limit. First, we derived 
the first order differential system for the geodesies in the black hole background. The result 
shows that the NUT charge n indeed influences the geodesies of particles. Then, combining 



- 20 - 



73 



.75 



a/2M=0.10 
a/2M=0.14 
a/2M=0.18 
a/2M=0.22 



7i \ 
76 \ 



74 



Figure 15. The first six caustics for tlic KTN black liole for different values of the spin a with fixed 
charge n = 0.08. 



with the lens equation and the null geodesies, we studied the equatorial lensing by the KTN 
black hole. The strong deflection limit coefficients o, h and the deflection angle a{u) were 
obtained. From the numerical calculation, we found that, for a fixed NUT charge n, a grows 
with the increasing of the spin a, while h decreases. With the values of the strong deflection 
limit coefficients a and 6, we got the numerical value of the deflection angle. In Fig. 6, it 
was plotted as a function of the spin a for different values of the NUT charge n. From it, 
the Kerr black hole and the KTN black hole were found to share the same property that 
the deflection angle a{u) monotonically increases with the spin a and the deflection angle of 
prograde photon is larger than the retrograde one. However, it was also found that, for a 
fixed value of spin a less than a negative value, the deflection angle a{u) for the Kerr black 
hole is smaller than that for the KTN black hole. And for a spin a greater than the value, 
a{u) of the Kerr black hole is larger than that of the KTN black hole. 

Considering that, even if the central ray of the light bundle is in the equatorial plane, 
most of the rays in the bundle will leave the equatorial plane, we investigated the quasi- 
equatorial lensing by the KTN black hole. We first got the precession of the orbit for small 
declinations for the KTN black hole spacetime, which is very different from the Schwarzschild 
black hole case. Then we studied the lensing at small declination. The strong deflection 
limit coefficients were reconsidered and their behaviors were found to be different from the 
equatorial ones. Under this approximation, the positions and magnification of the images 
were studied. Furthermore, the critical curves and caustic structure were obtained. The 
critical points are close to the optical axis for the positive spin a and farther for the negative 
spin a, and the nontrivial caustic structures acquire a finite extension. We also plotted the 
first six caustics for the KTN black hole for different spin a and NUT charge n, respectively. 
Compared all these results to those for the Schwarzschild black hole and the Kerr black hole, 
it is easy to find that there is a significant effect of the NUT charge n on the observables, 
i.e., the magnifying power /x/j, ratio fik+i/fik^ aiid the intersections Especially, for fixed 
value of spin a, the intersections 0™ vary from 0.1 ~ 2//arcsec for different value of NUT 
charge. Such an optical resolution is reachable by very long baseline interferometry (VLBI) 
projects, and advanced radio interferometry between space and Earth (ARISE), which have 
the angular resolution of 10^'^ arcsecond in the near infrared [68, 69]. Thus, measuring 6^ 
from astronomical observations, we are allowed to determine the value of the NUT charge n. 
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